Chaotic scattering on surfaces and collisional damping of collective modes 



Klaus Morawetz 1 ' 2,3 , Michael Vogt 4 , Uwe Fuhrmann 3 , Pavel Lipavsky 5 , Vaclav Spicka 5 
1 IPC, IN2P3-CNRS, ISMRA, 6 Bd du Marechal Juin, Caen, 14050, France 
2 GANIL, B.P. 5027, Bd Henri Becquerel, 14076 Caen Cedex 5, France 
3 Fachbereich Physik, University Rostock, D-18055 Rostock, Germany 
4 Institute of Philosophy, University Rostock, D-18055 Rostock, Germany 
5 Institute of Physics, Academy of Sciences, Cukrovarnickd 10, 16200 Praha 6, Czech Republic 



The damping of hot giant dipole resonances is investigated. The contribution of surface scattering is 
compared with the contribution from interparticle collisions. A unified response function is presented 
which includes surface damping as well as collisional damping. The surface damping enters the 
response via the Lyapunov exponent and the collisional damping via the relaxation time. The 
former is calculated for different shape deformations of quadrupole and octupole type. The surface 
as well as the collisional contribution each reproduce almost the experimental value, therefore we 
propose a proper weighting between both contributions related to their relative occurrence due to 
collision frequencies between particles and of particles with the surface. We find that for low and 
high temperatures the collisional contribution dominates whereas the surface damping is dominant 
around the temperatures y/3/2iv of the centroid energy. 



I. INTRODUCTION 

The damping mechanisms of collective motions in ex- 
cited nuclei are a topic of continuing debate Q . Mainly 
two lines of thought are pursued. In one line of thought 
it is assumed that collisions are the physical reason for 
damping only as developed via a Fermi liquid approach 
with bulk matter properties The other line of 

thought considers new features of the finite nucleus, such 
as surface oscillations and a level density with finite spac- 
ing. Partially the investigations are performed without 
inertia []l6|-p4| or by including inertia [§2 25-3(]]; note 
that inertia is absent in infinite matter. 

Both classes of models predict a comparable degree of 
damping necessary to reproduce the experimental data. 
Consequently, it is an open question which is the correct 
physical reason for damping. Of course, the correct de- 
scription has to assume a finite nucleus consisting of nu- 
cleons which are bound via the mean field, through which 
the nucleons undergo mutual collisions and where the sur- 
face is formed by the particles themselves. These features 
are usually included in Boltzmann-Uchling-Uhlcnbcck- 
(BUU) simulations |3l]-|33fl or in its nonlocal extensions 
p4,B5[. In full simulations, however, we will not gain a 
simple insight into the physical origins of the damping 
mechanism, in particular, how much is due to surface 
contributions and how much is due to collisional contri- 
butions. 

The aim of this article is to compare both pictures in 
the frame of linear response theory. Within the collision- 
free Vlasov equation the linear response of finite systems 
is well known |36| and allows one to calculate the strength 
function of finite nuclei. The damping, however, does not 
reproduce the experimental damping of giant resonances 



since collisions are absent. 

To include damping, we can take into account on 
the one hand the collisional damping in infinite matter 
and find a scaling for finite size effects in the sense of 
a Thomas-Fermi local density approximation. On the 
other hand we can consider the boundary of the finite 
nucleus as a fixed surface such that the nucleons sim- 
ply bounce off this wall like in a billiard-type situation. 
Provided we accept this simple picture, it is possible to 
compare the damping caused by this chaotic scattering 
off the wall with the damping from collisional contribu- 
tions in infinite matter. We will derive a response func- 
tion which includes the additional chaotic process and 
find a total damping rate T = T co n + r sur f similar to the 
Matthiessen rule in metals |37j but with an additional 
weight between the different processes according to the 
ratio of the corresponding collision frequencies. Moreover 
the hypotheses by Swiatecki |3£| is questioned that neg- 
ative curved surfaces of octupole deformations would in- 
duce an additional chaotic mechanism of dissipation and 
consequently octupole modes should be overdamped. We 
will find that the octupole deformation does not lead to 
any special enhancement of the Lyapunov exponent in 
comparison with the quadrupole mode. Moreover, any 
deformation will cause a contribution to the damping of 
collective modes. 

The outline of the paper is as follows. In the next chap- 
ter H the largest Lyapunov exponent for a nucle on i n a 
deformed nucleus is calculated. Then in chapter III we 
derive a unified response function which combines both 
the collisional damping and the contribution from surface 
scattering. In chapter IV both damping contributions are 



compared with the experimental values of hot isovector 
giant dipole resonances (IVGDR). A proper weighting 
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factor corresponding to the relative collision frequencies 
yields a unified picture which describes the data rather 
well. We find that the collisional contributions domi- 
nate for low temperatures and for very high tempera- 
tures while around the highest experimental achievable 
temperatures the surface effects are important. Chapter 
will summarize the results. 



II. LARGEST LYAPUNOV EXPONENT OF 
DEFORMED NUCLEI 

We will consider in the following only classical 3-D 
closed billiards and use the largest Lyapunov exponent 
as a relevant measure to characterize chaoticity. The 
surface of these billiards is chosen to resemble the sur- 
face deformation of a nucleus undergoing quadrupole or 
octupolc deformations. 

The Lyapunov exponent can be given as the deviations 
of the difference of the trajectories 
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as 



A = lim lim - In } J f 1 ^J, 

t^oo e^o t |rl(0) - r2(0)| 



(2) 



For A > the difference in phase space trajectories grows 
exponentially leading to chaotic behavior. 

We solve the Hamilton equations for one particle in 
finite nuclear matter with an infinite potential^] 



V(r,6) = V Q (r-R(6)) 



(3) 



modeling the deformation of a axial symmetric nuclear 
surface 



Rx(6) = R (1 + a 00 + a x P A (cos (9))) 



(4) 



with the nuclear radius Rq = 1.13^4 1 / 3 fm, and where 
A = 2 corresponds to the quadrupole and A = 3 to the 
octupole deformation p9[| . The coefficient aoo is adjusted 
to conserve the volume corresponding to incompressible 
nuclear matter. 



*The numerical implementation is performed with [Vb — > oo 
and a — > 0] 
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resembling an infinite step function. 
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FIG. 1. Overview of the geometry of modes 

The Lyapunov exponent is calculated by considering 
the time evolution of small deviations from a reference 
trajectory due to infinitesimal initial changes. We have 
used the Brandstaetter method resetting the deviations 
of the reference trajectory repeatedly after a certain time 
to the initial infinitesimal difference. This corresponds to 
an averaging and the largest mean Lyapunov exponent is 
obtained. 
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FIG. 2. The largest Lyapunov exponent of a spherical 
deformed billiard versus deformation parameter (^). The 
regions of a where the surface starts to become negatively 
curved (k < 0) via (jE)) are indicated by dotted lines at=o for 
octupole and quadrupole deformation, respectively. 

Figure ^| shows the Lyapunov exponent for different 
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octupole and quadrupole deformations according to (|J). 
The error bars are taken from averaging 50 runs of differ- 
ent initial conditions indicating 95% confidence level. We 
see that the quadrupole deformation leads to an imme- 
diate increase of the Lyapunov exponent while the Lya- 
punov exponent for the octupole deformation increases at 
larger deformation parameters and proceeds more slowly 
with increasing deformation. Let us note that the char- 
acteristic time scales from the Lyapunov exponent are 
similar to the damping rate of the IVGDR which will be 
presented in figure^. This motivates us to consider a uni- 
fied response from collisions and surface contributions. 

Following p8| a surface with negative curvature should 
induce a new chaotic mechanism analogous to the Sinai - 
billiard. We discuss therefore the curvature of a rotating 
body given by the curvature of the boundary curve R{9) 
of M\ via 



k = 



R 2 



2R 12 - RR h 



(R 2 + i?' 2 ) 3 / 2 



(5) 



where R' and R" denote the derivatives with respect 
to 9. In figure [2] we see that the octupole deformation 
shows parts with negative curvature at smaller deforma- 
tion parameter than the quadrupole deformation shows. 
However, the Lyapunov exponent of octupole deforma- 
tion remains always smaller. We therefore conclude that 
the octupole deformation compared with the quadrupole 
deformation leads to no significant enhancement of chao- 
tization. 



III. DAMPING OF COLLECTIVE OSCILLATIONS 

We would like to focus on a unified description of the 
linear response including collisional contributions and 
chaotic scattering with the surfaces. Therefore let us 
briefly sketch the response function formalism starting 
from appropriate kinetic equations. For infinite matter 
this procedure will result in the known Mermin response 
function which has been used to describe the IVGDR 
in symmetric JH| and asymmetric nuclear matter [ jl3| . 
Then we derive the response function for finite nuclei and 
show that in a local density approximation chaotic scat- 
tering from the surface can be incorporated. The result 
will be a Mermin-like response function in local density 
approximation where the largest Lyapunov exponent ap- 
pears as an imaginary shift in the frequency. 



A. Infinite matter response 

In order to consider the collisional damping we start 
from the kinetic equation for the quasiclassical distribu- 
tion function. For neutrons this kinetic equation reads 

P 

f n (p, r, t) H d r fn{p, r, t) - d r (V n + V cxt )d p f(p, r,t) = I 

m 



(6) 



with the self-consistent mean-field potential given by a 
schematic Skyrme type |40|,[fl| 



V n (r, t) = t Q { (f + y ) (n»(r, t) + n p (r, t)) 
n n (r, t) 



X 



+ t -f{(n n (r,t)+n p (r,t)) 2 -n 2 n (r,t)}. 



(7) 



The neutron and proton densities are n n — J T^tfis ft 



(2-ttR) 3 

and n p — J T^ffp fpi respectively. The parameters are 
x = 0.48, t = -983.4 MeV fm 3 and t 3 = 13106 MeV 
fm 6 . The kinetic equation and mean field for protons 
are obtained by interchanging the corresponding densi- 
ties. We restrict here to symmetric nuclear matter. Gen- 
eralizations to asymmetric nuclear matter can be found 
in p3| . First we sketch the results for the collision free 
(I = 0) or Vlasov equation and than we will take colli- 
sions into account. 



1. Collision free 

The collective effects without collisions (I = 0) can be 
obtained from the linearization of the Vlasov equation 
(^) with respect to an external potential V cx t- In the 
case of isovector oscillations, the density variation 5n = 
considered as the difference between proton 
and neutron densities. The linearization of the difference 
of kinetic equations for neutrons (^|) and protons leads to 



Sn = n (SV eKt + V Sn) 



with 1 41 



Vo = 2(14 - V p ) = -| (x + ~) - |no 



(8) 



(9) 



where no = J r^fp fo = 0.16 fm 3 is the nuclear satu- 
ration density and IIq is the Lindhard-function 



n|f f (g,u,) = 4 



dp' qd P 'fo{p') 

(2nh)3p'q ' 
oj — zO 



(10) 



Here we first consider the homogeneous equilibrium /o(p) 
or infinite matter. 

The response function n connects the induced density 
fluctuation Sn with the external potential via 



Sn = UV ei 



(11) 



and the response function follows from (|8j) for infinite 
matter 



n inf ( 9 ,^) 



(12) 
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The zeros ui = f2 + i7 of the dielectric function e(q, u) = 
in (|l2| ) determine the collective modes with the energy Q 
and the Landau damping 7 of the collective excitation. 
Within the Steinwedel- Jensen model for IVGDR @, the 
wave vector scales like q — where the nuclear radius 

is Ro = l.lSA^im. 



trajectories. It allows us to introduce the local density 
approximation as a first order memory effect in the tra- 
jectories. For this purpose we rewrite the Vlasov equa- 
tion (||) in a slightly different way. Introducing the La- 
grange picture by following the trajectory x(t),p(t) of a 
particle we linearize the Vlasov equation (|6|) according to 
f(x,p,t) = fo(x,p) +Sf(x,p,t) as 



2. Collisional model 

In order to take collisions into account as a further 
damping effect beyond Landau damping, we start from 
a kinetic equation analogous to (||) with an additional 
collisional term I\p, r, t] on the right hand side. In ]Tl| ] 
we have derived a collision integral in a non-Markovian 
relaxation time approximation 



dt 



Sf(x(t),p(t),t) = d p f d x(t ) SV 



(17) 



and get with one integration 
Sf(x,p,t) 

-2m f dt' [ dx'^S(x' - x{t')) df ^l; Xl) SV(x', t + if). 



dp 2 



-00 —00 



t _ 



(18) 



I(p,r,t) 



f(p, r, t)- f{p, r, t) 
-(t-t) 



dt 



(13) 



The density variation caused by the external potential is 
obtained as 



f f dp 3 

with the dynamical non-Markovian relaxation time 5n(x,uj) — —2ms / dx' / -, -^d v i f {p 2 ,x') 

j j (Zirhy 



T m (u)) 



1 

TB 



4 VttZV 



(14) 



The Markovian relaxation time is given by = 
^£aT 2 , where a is the averaged spin-isospin proton- 
neutron cross section. This collision integral holds for 
low temperatures compared to the Fermi energy. The 
non-Markovian relaxation time arises from the coupling 
of collective modes regarding two-particle scattering and 
consequently describes the effect of zero sound damp- 
ing. The local equilibrium distribution / is determined 
through the conservation of the local current. The lin- 
earization of the kinetic equation leads then to the ex- 
tended response function of Mermin 1 43 



1 - 



WT+J 



(1- 



(15) 



where the selfconsistency leads to the replacement of LlQ nf 
by II^ 1 in ([l2|). Generalizations to asymmetric nuclear 
matter can be found in Jl3[ . 

The energy and damping rates are determined by the 
zeros of the (Mermin) response function [jllj 



feo + i 7 ) = i-v Q (q)n™(q,n + i"f) = o. (16) 



Here the damping rate represents Landau and collisional 
damping. 



B. Finite matter response 

In the next step we will present the finite matter re- 
sponse function in terms of a memory integral over all 



dt'e- lt ' w 



7 S(x' -x(t'))SV(x',Lu) (19) 



dt 



where s denotes the spin-isospin degeneracy. Compar- 
ing this expression with the definition of the polarization 
function IIo, 



Sn(x, to) 



dx'n (x,x',Lj) 6V{x',uj), 



(20) 



we are able to identify the polarization of finite systems 
as 



Tlo(x,x' ,ui) = —2ms 



dp 3 



(2 7T ny 



-d p 2f (p 2 ,x') 



dt'e' 11 '" 



7 S(x'-x(t')). (21) 



dt 



Further simplifications are possible if we focus on the 
ground state /o(p 2 ) = ®(Pf ~P 2 )- The modulus integra- 
tion of momentum can be carried out and the Kirzhnitz- 
formula BJ45J for the polarization function appears 



n (x,x', u) = 

mspf(x) 
4tt 2 ^ 3 

mspf(x) 



A-k 2 Ti 




I "'■"•■si 

S(x' - as(0)) 



J dt'e- lt '" J ^S(x'-x(t')) 



(22) 
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This formula represents the ideal free part and a contri- 
bution which arises by the trajectories x(t) averaged over 
the direction at the present time n p pf = mx(0). In prin- 
ciple, the knowledge of the evolution of all trajectories is 
necessary to evaluate this formula. Avoiding the latter 
expense, we discuss two approximations which will give 
us an insight into the physical processes behind. First 
the most restrictive one shows how the local density ap- 
proximation emerges. As an extension we consider then 
the influence of chaotic scattering on a surface. 



1. Local density approximation 

The local density approximation appears from (J2S 
when two simplifications are performed. Introducing 
Wigner coordinates R = (x + x')/2, r = x — x' one has 
to assume 



1. gradient expansion 



p f (R+-)^p f (R) + 0(d R ) 



(23) 



2. expansion of the trajectories up to first order his- 
tory 

x' - x(t') S3 -r - t'x = -r- t'^n p + 0(t' 2 ). (24) 

TO 



C. Influence of chaotic scattering with surface on 
damping 

Now we want to assume an additional chaotic scat- 
tering which will be caused e.g. by the curved surface. 
In order to investigate this effect we add to the regular 
motion (E4J) a small irregular part Aa; 



x{t') n-r-t' 



Pf 



Ax. 



(28) 



This irregular part of the motion has the direction of the 
velocity, ri p , and lasts a time A t . During this time an 
exponential increase in phase-space occurs controlled by 
the largest Lyapunov exponent A. Therefore we assume 
[*' < 0] 



p f n p 



A t exp[-A(t' - A t )] + const. 



(29) 



Since we are looking for the upper bound of Lyapunov 
exponent we can take ( p9| ) at the maximum A t — — 1/A. 
Furthermore, in the case of vanishing Lyapunov expo- 
nent, the regular motion (|24|) should be recovered. This 
determines the constant. We obtain finally 



x(t') 



Pf, 



1 - cxp(-At') 



(30) 



With this ansatz one derives from (122) the result 



With these two assumptions we obtain from ([22|) after 
trivial integrations 



n™ A (q,R,u,) = 

{oo 
o 



(25) 



where C = m f^ ■ This can be further integrated with 
the help of 



d y S v ^- = arctan(ImC - iReC) -1 

y 



= 2iln 



i + C 
i-C 



+ 7T [sgn (1 + C) + sgn (1 - C)] 



ImC^O 

(26) 



leading to the standard Lindhard result for (|25|). We 
recognize the ground state result for infinite matter jlO| ) 
except that the Fermi momentum Pf{R) has to be un- 
derstood as a local quantity with respect to the density 



TL^ A (q,R,u J )=Il™ t ( q ,p f (R) 1 u;) 



(27) 



For extensions beyond the local density approximation 
see ilH. 



IL{q,R,u) = - 



mspf(R) 



J y V w 





iu j A — 1 



(31) 

which for A — > resembles exactly (^5|). The further 
integration could be given in terms of hypergeometric 
functions but this is omitted here. 



With formula (31) we have derived the main result of 
this section, i.e. a polarization function due to many 
particle effects including the influence of an additional 
chaotic process characterized by the Lyapunov exponent 
A. 

For small values of -^A = C— which corresponds to 
relative small Lyapunov exponents, we can use lim (1 + 

x — >oo 

a/x) x = exp(a) in the integral of (|3l| ) and the final in- 
tegration can be performed with the result of (^?j) and a 
complex shift in the frequency 



ILT (q,R,u) = n-' (q, p f OR), to + iX) 



(32) 



Solving the dispersion relation we obtain in this way 
the Matthiessen rule which states that different damping 
mechanisms are additive in the total damping. However, 
the collisional damping of infinite matter, 27, and the 
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largest Lyapunov exponent, A, are not yet the appropri- 
ate values of finite matter. Instead, these two damping 
mechanisms have to be added by a proper relative weight 
which will take into account the relative occurrence of 
the processes, the collisions with other particles and the 
collisions with the surface. 



IV. COMPARISON OF DAMPING 
MECHANISMS 

A. Connection between surface curvature and 
temperature 

The damping of excited nuclei can be understood by 
two different mechanisms. Besides the collisional damp- 
ing we also have to consider the shape fluctuations. In the 
absence of inertia which will be assumed in the following, 
the driving force for shape fluctuations is the tempera- 
ture. Therefore we link the surface deformation a of (Q) 
to the temperature within a statistical model. We use as 
a measure for the mean deformation F 



(a) 



j da\a\exp(—EB(a)/T) 
J da exp(-E B (a)/T) 



(33) 



where the surface dependent energy Es(a) is given by 
the Bethe-Weizsiicker formula ^ 



Eb{ol) = —ai + 



a-2 



0*&_ 

A l/3 + A i/3 



a4<5 + a^A 



2/3 



S(a) 
5(0) 
(36) 



with the volume energy a\ = 15.68 MeV, Coulomb en- 
ergy a>3 = 0.717 MeV, the symmetry energy 04 = 28.1 
MeV and the surface energy 05 = 18.56 MeV. 

The surface of rotational symmetric nuclei (^) is given 

by 



2tt 



dxR x {x)^R x {xf + (1 - x 2 )R' x { X y (37) 



^For small deviations we found identical Lyapunov expo- 
nents for prolate a > and oblate a < deformations and 
therefore we do not distinguish the sign of a. 

* Please remind that in principle the Coulomb energy ch ang es 
with small deformation as well according to the factor Wm 

1 ™ 

while the surface term changes as 

1 + (A- l)(A + 2)/2/(2A + l)ai (35) 

Only the correction ( [35| ) is considered since (Q) would lead 
to corrections of around 0.3% and are neglected here. 



from which we obtain for the quadrupole, S2, and oc- 
tupole, S3, deformations (B5|) 



(38) 



This represents the lowest order expansion in a, however 
the next term gives already corrections in fractions of per- 
cents for the highest deformations considered here. By 
this way, the statistical model ([53]) leads to a connection 
between temperature and deformation as 



T : 



.2/3 



(39) 



where the constant c is given by the coefficient of a 2 in 
Eq. (38) for the corresponding quadrupole or octupole 
deformation. 



B. Finite matter scaling for collisional damping 



According to the local density approximation ( |27| ) we 
want to investigate the collisional damping in finite mat- 
ter. Therefore we have to replace all densities in the 
dispersion relation ( |l5| ) by the local density which are 
parameterized by a Woods-Saxon potential 



n(r) 



n 



exp[(r-i? )/0.545fm] + l' 



(40) 



Solving the dispersion relation ( fl6|) w ithin the local den- 
sity approximation according to (p7[) we obtain a spatial 
dependent damping rate 7(7"). 

In the following we derive a finite size scaling factor, 
£a , of the damping 



<7(0) 
7 



(41) 



which approximates the local (r-dependent) dispersion 
relation by an averaged one 



(...) 



ft 3 , 

-^o Jo 



R 



d r 



(42) 



assuming radial symmetry. The factor ( pd"| ) has to be 
applied to the infinite matter bulk collisional damping 
and is found &n = 0.9199 and £ Pb = 0.9250 for 120 Sn 
and 208 Pb respectively. Therefore, the bulk matter value 
of damping is diminished by finite size effects. 



C. Comparison of collisional damping and surface 
damping 



Using Eq. fl39| ) we can translate the Lyapunov expo- 
nent A calculated as a function of deformation in figure 
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into a function of the temperature. In figure^ the contri- 
bution to the damping of IVGDR for 120 Sn (circles) and 
208 Pb (squares) is presented for different shape deforma- 
tions versus temperature. If we add quadrupole and oc- 
tupole deformations we come up with a damping curve 
very similar to the paper of Ormand [^2| . The damp- 
ing starts at zero and increases rapidly with increasing 
temperature. We see that the main contribution comes 
from the quadrupole deformation while the octupole de- 
formation is only sizeable at higher temperature. Let us 
note that the qualitative difference between Sn and Pb 
is reproduced by surface scattering as well as collisional 
damping. 




temperature T [MeV] 

FIG. 3. The collisional damping scaled with finite size ac- 
cording to ( pl| ) is compared with the damping according to 
the chaotic scattering from the surface of quadrupole and oc- 

120 208 

tupole deformed shapes for Sn and Pb from figure El 

The collisional contribution 27 scaled to finite sizes via 
(fy]) are plotted as well (Sn: solid line, Pb: dashed line). 
We recognize that both contributions by itself, collisional 
as well as surface scattering, account almost for the same 
amount required by the experimental values, see figure^]. 
A proper relative weight between both processes is there- 
fore necessary which will be introduced in the following. 

So far we have not considered that only particles close 
to the surface can appreciably contribute to the sur- 
face chaotization, while particles deep inside the nuclei 
are screened out of this process. Consequently we con- 
sider the corresponding collision frequencies as the mea- 
sure to compare surface collisions with interparticle colli- 
sions. The collision frequency between particles is given 
by 1/tm of (14). The collision frequency of particles 
with the deformed surface beyond a sphere, v S nil , is given 
by the product of the density with the surface increase 
S(a) - S(0) = coPA-kRq according to Eq. ( p8| ) and with 
the mean velocity in radial direction v r = 3/8vf- The 
result is 

^surf = 1.5Tn Q v F rl/a 5 (43) 
where we have used Eq. (p9) to replace a. We see that 



the frequency (|43|) is independent of the size of the nu- 
cleus and linearly dependent on the temperature. 

We use the ratio of these two frequencies to weight 
properly the two damping mechanisms, the surface colli- 
sional, A, and interparticle collisional, < "f(r) >, contri- 
butions. Consequently the full - width - half - maximum 
reads (FWHM) 



Tfwhm = 2 (C < j(r) > +(1 - C) A) 

= r co n -(- r sur f. 



With the help of (|14j) and 
given by 

C(T) = - 



(44) 



fj the weighting factor £ is 



T M {T) 



(45) 



One sees that for zero and high temperatures £ = 1 and 
due to Eq. ([l4|) only the collisional contributions matter. 
Since u sm { is linear in the temperature and 1/r depends 
quadratically on the temperature, the weighting factor £ 
has a minimum at temperatures around T c = and 
the surface contributions become important. In the case 
of the IVGDR this corresponds to a temperature of T ks 
3.7MeV, which is the upper limit of current experimental 
achievable temperatures. Therefore we can state that 
at low and high temperatures the collisional damping is 
dominant while for temperatures around T c the surface 
contribution becomes significant. 

In figure w e compare the effective damping accord- 
ing to Eq. (pd|) with the experimental data. We find a 
reasonable quantitative agreement. 




temperature T [MeV] 

FIG. 4. The effective damping consisting of collisional and 
surface damping to geth er with the experimental data (filled 
symbols) from Ref. |TJ] (Sn) and from Ref. ( pb )- 

While the collisional damping by itself can reproduce 
the increase at small temperatures well enough it fails 
to reproduce the sharp rise at higher temperatures, see 
figure ^. The surface damping by itself fails at zero tem- 
perature and leads to an unphysical shape. Only the 
connection of both effects, collisional and surface damp- 
ing seems to be able to reproduce the data correctly. 
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At this point we like to stress that we have adopted 
a simplified model of finite Fermi liquids and introduced 
a phenomenological weighting between surface and in- 
terparticle collision contributions proportional to their 
relative occurrence. The discussion of damping due to fi- 
nite size effects within shell model considerations can be 
found in literature, e.g. the review where this problem 
has been addressed from the point of doorway states. 

V. SUMMARY 

We have considered the influence of collisional damp- 
ing as well as surface scattering for the damping of hot 
isovector giant dipole resonances (IVGDR) . We find that 
the relative importance of both effects is dependent on 
the temperature scale. While for low and high tempera- 
tures the collisional contribution dominates, the surface 
scattering is the one of significant importance at temper- 
atures around 28% of the centroid energy. 

The surface contribution has been calculated from the 
Lyapunov exponent of surface scattering for different de- 
viations of the nucleus shape. The deviations are related 
to the temperature by a statistical model. 

The influence of both the surface and collisional damp- 
ing is described by a generalized response function for fi- 
nite nuclear matter which takes into account the chaotic 
processes of scattering with the surface. We derive by 
this way a response function similar to the Lindhard re- 
sponse in local density approximation which is now mod- 
ified by the Lyapunov exponent of the surface scattering. 
In an approximative way we show that the total damping 
is described by the sum of the two contributions to the 
damping in analogy to the Matthiessen rule. 

Comparing the collision frequencies for one particle 
with other particles and for one particle with the surface, 
we derive a proper relative weight for both processes, the 
collisional and the surface damping. The resulting effec- 
tive damping reproduces the experimental data rather 
well in the whole accessible range. 
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